Direct detection of the inflationary gravitational wave background 
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Inflation generically predicts a stochastic background of gravitational waves over a broad range 
of frequencies, from those accessible with cosmic microwave background (CMB) measurements, to 
those accessible directly with gravitational-wave detectors, like NASA's Big-Bang Observer (BBO) 
or Japan's Deci-Hertz Interferometer Gravitational- wave Observer (DECIGO), both currently under 
study. Here we investigate the detectability of the inflationary gravitational-wave background at 
BBO/DECIGO frequencies. To do so, we survey a range of slow- roll inflationary models consistent 
with constraints from the CMB and large-scale structure (LSS). We go beyond the usual assumption 
of power-law power spectra, which may break down given the 16 orders of magnitude in frequency 
between the CMB and direct detection, and solve instead the inflationary dynamics for four classes of 
inflaton potentials. Direct detection is possible in a variety of inflationary models, although probably 
not in any in which the gravitational-wave signal does not appear in the CMB polarization. However, 
direct detection by BBO/DECIGO can help discriminate between inflationary models that have the 
same slow-roll parameters at CMB/LSS scales. 

PACS numbers: 98.80.Bp,98.80.Cq,04.30.Db,04.80.Nn 
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I. INTRODUCTION 



Long a subject of theoretical speculation, inflation 
IE 0] has now, with the advent of precise cosmic mi- 
crowave back grou nd (CMB) measurements 0, IE IE IE 
IE IE EE EH IH] become an empirical science. The 
concordance of the measurements with the inflation- 
ary predictions of a flat Universe and a nearly scale- 
invarian t sp ectrum of primordial density perturbations 
[l3L IhH Il5l fl(| is at least suggestive and warrants fur- 
ther tests of inflation. Among the predictions of in- 
flation yet to be tested is a stochastic gravitational- 
wave background with a nearly scale-invariant spec- 
trum IH HE EE IE El EE O. Detection of the 
CMB polarization pattern induced by inflationary gravi- 
tational waves of wavelengths comparable to the horizon 
has become a goal of next-generation CMB experiments 
ElEEEEElllil. And now, direct detection of the in- 
flationary gravitational wave background (IGWB) with 
future spaced-based gravitational- wave detectors at deci- 
Hertz frequencies has become the subject of serious study 

mm. 



Detection of a gravitational-wave background, at ei- 
ther CMB or direct-detection frequencies, would consti- 
tute a "smoking gun" for inflation. Moreover, since the 
amplitude of the IGWB is determined by the energy scale 
of inflation at the time that the relevant distance scale ex- 
ited the horizon during inflation, detection would provide 
important information about the new ultra-high-energy 
physics responsible for inflation [3ll l32| . Since the fre- 
quencies probed by the CMB and by direct detection are 
separated by 16 orders of magnitude, the combination of 
both provides a large lever arm with which the shape of 
the inflaton potential can be constrained. 

In this paper, we survey a range of inflationary mod- 
els to investigate the detectability of the IGWB with 



satellite experiments, like NASA's Big Bang Observer 
(BBO) [23 and the Japanese Deci-Hertz Interferome- 
ter Gravitational-wave Observatory (DECIGO) [3(|, cur- 
rently under study. We restrict our attention to slow-roll 
inflation models that are consistent with measurements 
from the CMB and large-scale structure. We show how 
measurements of the IGWB amplitude at both CMB and 
direct-detection scales can be used to constrain the infla- 
tionary parameter space. 

Previous work [33l l34l EE] on direct detection of the 
IGWB has taken the gravitational-wave spectrum to be 
a pure power law, considered chaotic inflation [2E E3 
or the IGWB due to a broken scale invariant potential 
|38|. In this paper, we consider a wider range of in- 
flationary models (in the spirit of Refs. |3E E3)> an d 
we solve the inflationary dynamics to go beyond the as- 
sumption of power-law power spectra. With this more 
accurate analysis, we find that for the forms of the infla- 
ton potential considered here the direct detection of the 
IGWB can break degeneracies between distinct inflation- 
ary models that produce the same slow-roll parameters 
at CMB/large-scale-structure scales for broken scale in- 
variant potentials. 

This paper is organized as follows. In Section II, 
we review the basics and relevant ingredients of infla- 
tion as well as constraints from the CMB and large- 
scale structure (LSS) to the inflationary observables. We 
also discuss the transfer function that relates the current 
gravitational-wave amplitude to its primordial value. In 
Section III, we discuss the sensitivities of planned space- 
based gravitational- wave observatories to the IGWB. In 
Section IV, we show for four different families of slow-roll 
inflation models the IGWB amplitude and spectral in- 
dex rit at BBO/DECIGO scales that are allowed by cur- 
rent CMB/LSS constraints to the n s -r parameter space 
(where n s is the scalar spectral index and r the tensor- 



2 



to-scalar ratio). Section V compares the results of our 
calculations with those obtained by extrapolating the 
power-law power spectra inferred from CMB/LSS mea- 
surements to BBO/DECIGO scales. In section VI we 
discuss a family of slow-roll inflation models where the 
observational signature in CMB/LSS studies for two dif- 
ferent models is the same but differs for the direct detec- 
tion of the IGWB at BBO/DECIGO scales. In Section 
VII, we summarize and make some concluding remarks. 



II. INFLATIONARY DYNAMICS AND 
PERTURBATIONS 

A. Homogeneous Evolution 



the usual slow-roll parameters, 
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which are required to be small compared with unity for 
the slow-roll approximation to be valid. Toward the end 
of inflation, e grows, and inflation ends when e ~ 1. This 
statement can be made precise by the use of "Hubble 
slow-roll" parameters |4lj . 



Inflation occurs when the cosmological expansion ac- 
celerates; i.e., when a > 0, where a{t) is the scale factor, 
and the overdot denotes a derivative with respect to time 
t. The evolution of the scale factor is determined by the 
Friedmann equation, 



H z = [ - 
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the continuity equation, p + 3H(p + P) = 0, and an equa- 
tion of state P(p), where H is the Hubble parameter, p 
is the total energy density, P is the pressure, Mpi is the 
Planck mass, and K is a constant related to the 3-space 
curvature. From Eq. and the continuity equation fol- 
lows the "acceleration" equation, 
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For an equation of state of the form P = wp, where w is 
a constant, inflation occurs when w < —1/3. 

Consider now a spatially homogeneous scalar field 0, 
the "inflaton" . It has an energy density and pressure, 



1 



P = 



P = 



v(4>), 



(3) 
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from which it follows that inflation occurs if V{<j>) > <j) 2 . 

The equation of motion for the inflaton is given by 
+ 'iH(f> + V = 0, where the prime denotes differenti- 
ation with respect to (p. We assume that inflation has 
been proceeding for a long time before any observable 
scales have exited the horizon, and so for our purposes, 
the energy density is dominated by the inflaton during 
inflation, the curvature term, K/a 2 , is negligible as com- 
pared to the inflaton energy density, and the evolution 
of the inflaton has been attracted to the slow-roll regime 
(e.g., Ref. 01 )■ ^ so ' the evolution of the inflaton and 
the scale factor are uniquely determined by V(</>). Within 
the slow-roll approximation, the evolution is described by 



B. Perturbations 

To leading order in the slow-roll approximation, the 
amplitudes of the power spectra for density perturbations 
(scalar "s" metric perturbations) and gravitational waves 
(tensor "t" metric perturbations) can be written (e.g., 
Refs. HHH) 
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as a function of wavenumber k, where V and V are eval- 
uated when the relevant scale exits the horizon during 
inflation. The power spectra can be expanded in power 
laws, 
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nt + (at/2) ln(fc/fco) 
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where fco is a pivot wavenumber at which the spectral 
parameters (e.g., Ref. [H]). 



n s (fc) ~ I- 66 + 27?, 

n t {k) ~ -2c, 

a s (k) ~ 16677 - 24e 2 - 2£, 

a t (k) ~ 4e7y - 8e 2 , 
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are to be evaluated. To a first approximation, the power 
spectra are power laws with power-law indices n s and n t , 
although these indices may "run" slightly with k, with a 
running parameterized by a s and a t |45j . Finally, the 
tensor-to-scalar ratio is 



P(k) 
P s (k) 



16e. 



(16) 
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In this paper, we will generally evaluate P s , n s , and a s at 
the distance scales of the CMB and large-scale structure 
(LSS), where they are measured or constrained. In the 
Figures below, the tensor spectral index nt will be eval- 
uated at the distance scale relevant for direct detection 
of gravitational waves. 



C. Number of e-Foldings 

The number of e-foldings of expansion between the 
time, determined by k — a^Hk, when a comoving dis- 
tance scale labeled by k exited the horizon during in- 
flation, and the end of inflation is N{k) = ln(a en d/afc), 
where a en d is the scale factor at the end of inflation. This 
must be (e.g., Ref. |4f|), 



N{k) 
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where is the inflaton potential when the comoving 
scale k crossed the inflationary horizon; V en d is the value 
of the potential at the end of inflation; and p r h ~ T^ h is 
the energy density of the universe once radiation domi- 
nation begins. Here we have assumed that the Universe 
was matter dominated after inflation and before reheat- 
ing and that the transition between radiation and matter 
domination is instantaneous. In terms of the inflaton po- 
tential, the number of e-foldings between two field values, 
4>i and </>/, is 



N(4n,<f>f) 
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where we have supposed the potential increases as the 
field increases so that the field rolls towards the origin. 
Furthermore, if the potential determines a field value at 
which inflation ends, we can combine this equation with 
Eq. I|17|) and p r h to identify the field value when the cur- 
rent Hubble volume exited the inflationary horizon. 

If we use 10 16 GeV for all the densities in Eq. l(T7|l . we 
require 62 e-folds of inflation between the time the cur- 
rent horizon distance exited the horizon during inflation 
and the end of inflation. The strength of the IGWB is 
proportional to the inflaton-potential height [cf., Eq. ©], 
and as we will see, detectability requires V > 10 15 GeV. 
We will also see that V cn d is never much smaller than 14 . 
Thus, the only ratio in Eq. (|17|l that might be large is the 
last. Conservatively, the reheat temperature must be > 
1 MeV to preserve the successes of big-bang nucleosyn- 
thesis, implying a lower limit of N(aoHo) > 47. This 
lower limit is significant, as the ratio of gravitational- 
wave frequencies probed by the CMB/LSS (correspond- 
ing to k ~ 0.05 Mpc -1 ) and direct detection is ~ e 35 . 
Therefore, if inflation results in an IGWB in the ballpark 
of detectability by the CMB, then inflation will last long 
enough to ensure the production of gravitational waves 



on BBO/DECIGO scales (although it does not necessar- 
ily guarantee a detectable amplitude). 



D. Constraints to Inflationary Observables 

We would like to survey only those inflationary mod- 
els that are consistent with current data. Measure- 
ments of the "inflationary observables" — i.e., the scalar 
and tensor power-spectrum amplitudes, spectral indices, 
and running — come from CMB measurements that probe 
scales from the current Hubble distance (~ 10 4 Mpc) 
to ~ 10 Mpc scales, galaxy surveys that constrain the 
matter power spectrum from 600 Mpc down to 20 Mpc, 
and from the Lyman-a forest, which probes down to ~ 
1 Mpc 0! Oil • Constraints to the classical cosmologi- 
cal parameters (e.g., the Hubble parameter, the deceler- 
ation parameter, the baryon density, the matter density) 
from other measurements help limit the range of plausi- 
ble values for the inflationary observables that come from 
CMB/LSS measurements. 
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FIG. 1: Regions in the n s -r parameter_space consistent with 
the CMB-only (orange/medium gray) |49l| . CMB plus galaxy 
surveys (red/dark gray), and CMB plus galaxy surveys plus 
Lyman-alpha- forest constraints (green/light gray) 15011 . Here, 
r is the tensor-to-scalar ratio, and n s is the scalar spectral 
index at CMB scales. We plot on top of these regions the 
parameter space occupied by the four models of inflation we 
consider: power-law (solid black line), chaotic (dotted red), 
symmetry-breaking (dashed-dot purple) and hybrid (short- 
dashed dark-red). The right axis shows the energy scale 
[V(k c )] 1/4 of inflaton. 

The precise constraints to the inflationary observables 
depend in detail on the combination of observational data 
sets. In our discussions, we simply take as conservative 
ranges P s (k ) = (2.45±0.23) x 10~ 9 , n s = 1.0±0.1, and 
\a s \ < 0.04 at a pivot wavenumber ko = 0.05 Mpc -1 |Hc| . 
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We note that the CMB-only constraints correspond to a 
pivot wavenumber, k = 0.002 Mpc -1 49]. The errors we 
quote are not to be interpreted as statistical errors bars, 
except in the case of the value for P s (ko), where the er- 
ror quoted is Ict; rather, they simply indicate a range of 
parameters that are in reasonable concordance with mea- 
surements 0, El an d the range of values that 
we use here. In our discussion, we take a conservative 
upper limit r < 1 to the tensor-to-scalar ratio. The nu- 
merical results we show in Fig.[3]use slow-roll parameters 
consistent with the regions in the n s -r parameter space, 
shown in Fig. ^ taken from analyses of CMB data alone, 
CMB plus galaxy surve ys, and CMB, galaxies, and the 
Lyman-alpha forest |49ll5(]j]. 



E. Gravitational- Wave Transfer Function 

The gravitational-wave power spectrum Pt(k) provides 
the variance (|/ifc| 2 ) of the gravitational-wave amplitude 
hk as that mode enters the horizon. Once the wavelength 
is smaller than the horizon, the gravitational wave begins 
to oscillate, and its energy density pk ~ k 2 h\ redshifts 
oc a~ 4 like radiation. It follows that the gravitational- 
wave amplitude today is hk(to) = hk{tk){dk/ao), where 
to is the time today, t k is the time of horizon entry, and 
a k = a(tfc) |103| . During radiation domination (RD), 
H oc a -2 , so that k oc and during matter domina- 

tion (MD), H oc a~ 3 / 2 , so that k oc 1/a^ 2 . From these 
relations, we find that the value of hk today scales with 
k as 

hk oc (horizon entry during RD), (19) 
hk oc k~ 2 (horizon entry during MD). (20) 

Calculations (e.g., Refs. 01 tne transfer func- 

tion intended for CMB predictions evolve the wave equa- 
tion more carefully through matter-radiation equality, 
but the direct-detection frequencies are so high that the 
scalings we have used here are fairly precise. The sen- 
sitivities of BBO/DECIGO will peak near a frequency 
0.1 Hz, or wavenumber k — 6.47 x 10 13 Mpc _1 . Matter- 
radiation equality corresponds to fc cq « 0.05/i 2 Mpc -1 . 
Therefore, the primordial gravitational waves observed 
by the planned gravitational-wave observatories entered 
the horizon well before matter/radiation equality. In 
fact, the modes that entered the horizon during big- 
bang nucleosynthesis (at T ~ 1 MeV) have frequen- 
cies ~ 10~ n Hz. Therefore, the gravitational waves 
probed by BBO/DECIGO are ~ 10 orders of magnitude 
smaller than those associated with big-bang nucleosyn- 
thesis (BBN) and must have entered the horizon at tem- 
peratures T ~ 10 7 GeV. 

These high temperatures imply a small correction to 
previous calculations, which assumed T oc a -1 , of the 
transfer function due to the fact that it is actually 
g*s{T)a T 3 that remains constant, where g*s(T) is the 
effective number of relativistic degrees of freedom con- 



tributing to the entropy density. Recalling that the Hub- 
ble parameter is H ~ 1.66 g*T 2 /Mj>i during radiation 
domination (where g* is the effective number of rela- 
tivistic degrees of freedom contributing to the energy 
density; at these temperatures, ~ g*s)> the condi- 
tion for horizon entry for a physical wavenumber today 
(aofc = akHk), becomes 



i .(.6 ff y 2 (T fe )(ToT fe /M P1 ) ( ^£f0 = 



(21) 



Today, g*s(Tb) = 3.91 for photons and three species of 
massless neutrinos. We then have, 



3.8 x 10 6 GeV 



g*(T k ) 



100 



1/6 



6.5 x 10 13 Mpc _1 



(22) 

Taking only standard-model particles, g*(T) w 100 (or 
roughly double if there is low-energy supersymmetry) 
and should be roughly independent of temperature. We 
thus find that for BBO/DECIGO scales, 



T(k) = -=2.1x 10" 
a Q 



6.5 x 10 13 Mpc 



-1/6 

(23) 

where f/ioo = 5*(T\.)/100. With this transfer function, 
the root- variance of the IGWB today is 



<IM 2 ) 1/2 =Pt{k) 1/2 T{k). 



(24) 



Free streaming of neutrinos, which occurs after neutri- 
nos decouple at a temperature ~ 1 MeV, contributes an 
anisotropic stress |55l l56| . However, the resulting damp- 
ing will be negligible for modes that enter the horizon so 
much earlier than neutrino free-streaming. 

If we define the logarithmic GW contribution to the 
critical density, 



^gw(A) = 



1 dp GW 
p r din k ' 



then \51 



ft G w{k)h 2 



6H 2 



■(\hk\ 



A G wPt(k), 



(25) 



(26) 



where H = 100 h knis -1 Mpc" 1 , and A GW = 2.74 x 
10 Sioo ■ l n slow-roll inflation, n t < 0, so if we take 
P s (ko) ~ 2.45 x 10 -9 and r < 1 (corresponding to V 1 ^ 4 < 
3.36 x 10 16 GeV), we find an upper limit Q GW h 2 < 6.71 x 



III. DIRECT-DETECTION THRESHOLDS 

Since the detection of a stochastic background of grav- 
itational waves has to be separated from the effect of 
sources of noise intrinsic to the detector, the sensitivity 
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to a stochastic background is different than the sensitiv- 
ity to a non-stochastic source. The total stochastic signal 
in a given detector can be written as a sum of a stochastic 
signal plus a stochastic noise, s(t) = h n (t) + h(t). Taking 
the Fourier transform of the total signal and considering 
the spectral density of the noise and of the stochastic sig- 
nal, we find that in order to have a signal-to-noise greater 
than unity, 




where F is a filling factor that accounts for the fact that a 
primordial stochastic background will be isotropic on the 
sky, but the detector will only be sensitive to a fraction 
of the sky, while S n (f) is the spectral density associated 
with the detector noise. 

For omni-directional interferometers, such as the Laser 
Interferometer Space Antenna (LISA), F = 2/5. There is 
a great improvement in sensitivity when the signal from 
two independent detectors can be combined th roug h a 
correlation analysis between the two detectors [E3TIEE 
IHqL IfSflj . Such a correlation increases the sensitivity to a 
stochastic background such that, 

(IM/)| 2 > 1/2 > (28) 




where A/ is the bandwidth over which the signals can 
be correlated and AT is the integration time. For a 
correlation analysis, the increase in sensitivity is under 
the assumption that the detector noises are indepen- 
dent between the two detectors, while the only corre- 
lation expected is due to stochastic signals such as in- 
flation. For the single detector, the minimum observ- 
able strain is independent of the integration time, while 
for a correlation analysis, long-term observations are es- 
sential. While LISA will not allow an opportunity for 
such a correlation analysis, some mission concept stud- 
ies for NASA's Big Bang Observer (BBO) and Japan's 
Deci-Hertz Interferometer Gravitational Wave Observa- 
tory (DECIGO) consider two (or more) systems such that 
the improvement related to the correlation analysis can 
be exploited. The design for LISA currently places the 
sensitivity at approximately f2cw^ 2 ~ 10~ n . Current 
designs for BBO place the sensitivity of a single detector 
at ficw^ 2 ~ 10 -13 and the sensitivity of a correlated 
extension at f^cw^ 2 ~ 10 -17 . Finally, the ultimate goal 
for DECIGO is a sensitivity to fl GW h 2 ~ 1CT 20 H3, cor- 
responding to V 1 / 4 ~ 1.15 x 10 15 GeV (r ~ 10" 6 ). 

Besides a sensitivity to a stochastic background, one 
must also be concerned about sources of a stochas- 
tic background, other than inflation. One such source 
is the background in extragalactic supernovae [6l| . 
Such sources have the potential to wash out any sig- 
nal that would otherwise be observed from a primor- 
dial source, but the characterization of the amplitude 



and frequency dependence of these sources is still un- 
certain. Other sources of cosmological gravitational- 
wave backgrounds are white-dwarf/whitc-dwarf binaries 
|62|. neutron-star/neutron-star binaries |63| and neutron- 
star/white-dwarf binaries |r3^|. 

Fig. [21 shows the sensitivities to a stochastic 
gravitational- wave background as a function of frequency 
/ for a variety of gravitational- wave detectors. Also 
shown are various current limits (as solid curves) as 
well as a variety of projected direct and indirect sensi- 
tivities (dashed curves), and scale- invariant spectra pa- 
rameterized by an energy scale V 1 ! 4 of inflation (dot- 
ted curves). We also show limits from current CMB ex- 
periments as well as the sensitivities expected for future 
CMB-polarization experiments currently under study. 

IV. BBO/DECIGO AMPLITUDES 

In this Section, we calculate the gravitational- wave am- 
plitude at BBO/DECIGO scales for several families of 
slow- roll inflation models consistent with CMB/LSS con- 
straints. 

Measurements of the scalar amplitude P s and spectral 
index n s at CMB/LSS scales, as well as upper limits to 
the tensor contribution r to the CMB and to the running 
a s of the spectral index, constrain the inflaton potential 
and its derivatives at the field value 4> c that corresponds 
to the time at which CMB/LSS scales k c exited the hori- 
zon. To be precise, we use k c = 0.05 Mpc -1 . In this 
work, we take as the nominal BBO/DECIGO frequency 
/ = 0.1 Hz, corresponding to k = 6.47 x 10 13 Mpc _1 
(and we note that flGw(k)h 2 ~ constant for n t ~ 0, so 
our results will not depend too sensitively on the precise 
value of / we use). CMB/LSS and BBO/DECIGO scales 
are therefore separated by AN = ln(6 x 10 13 /0.05) ~ 35 
e-folds of inflation 104]. Eq. i|18|) can then be used to 
find the field value (j) g at the time that BBO/DECIGO 
scales exited the horizon. 



A. Power-Law Inflation 

In power-law inflation, the inflaton potential takes the 
form, 

V{4>) = V e- p ^ Mpi . (29) 

Power-law inflation is so called because the scale factor 
is a power law a(t) cx t 16 */P , and the Hubble param- 
eter also evolves as a power of time t. The resulting 
scalar and tensor power spectra are then pure power laws, 
with no running of the spectral indices. The parameter 
e = p 2 /(16n) always, so that inflation must be ended ar- 
tificially at some </> e nd- Although the potential has only 
two free parameters (Vo and p), there is an additional 
free parameter, namely, the value of <p c , which we are 
free to choose in this particular family of models. This 
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FIG. 2: Current limits and projected sensitivities to a stochastic gravitational-wave background versus the gravitational-wave 
frequency. The solid curves all indicate current upper limits, while the dashed curves indicate projected sensitivities. The LISA 
curve is from Ref. |6Sl| and BBO correlated from Ref. |6l| . The BBO sensitivity is estimated by increasing the BBO-correlated 
curve by 4 orders of magnitude [see Eq. 11291 1. The BBN constraint results from the limit to the number of relativistic degrees of 
freedom at big-bang nucleosynthesis (e.g., Ref. pal): the "M/R" constraint is from CMB/LSS constraints to matter-radiation 
equality |67|: the "z. var" curve is from Ref. [68l|: and the quasar-astrometry limit from Refs. |68l. ITOB . We note that the BBN 
and "M/R" constraints assume a scale invariant gravitational- wave background that extends ~ 60 e- folds below the current 
Hubble horizon. LIGO sensitivities, taken from the LIGO Scientific Collaboration White Paper on Detector Research and 
Development [Till are given in terms of a correlated analysis between the Hanford, WA and Livingston, LA sites [see Eq. 129H 1. 
The run 1 LIGO limit ("SI LIGO") is from Ref. El and the run 3 LIGO limit ("S3 LIGO") is from Ref. [73|. Also shown 
are millisecond-pulsar timing constraints (current |74l |75| and sensitivities projected for the Square- Kilometer Array |76j|L 
Curves corresponding to scale-invariant (i.e., n t = 0) gravitational- wave backgrounds are shown (dotted curves), labeled by the 
associated inflationary energy scales at CMB/LSS scales (but keep in mind that slow-roll inflation generically predicts n t < 0, 
less power on small scales). The CMB/LSS currently constrains this value to be below 3.36 x 10 16 GeV at CMB/LSS scales. 
Future CMB measurements may be able to reach energy scales near 10 15 GeV Ellzilzlli^- 



model has also r\ = p 2 /8ir, so n s = 1 — p 2 /8ir = 1 — 2e, 
and for n s > 0.9 we find a constraint e < 0.05. The con- 
straint r = 16e < 1 is comparable or a bit weaker. Since 
n s and r depend in this model only on the parameter 



p, these models occupy a curve in the n s -r parameters 
space, which is indicated by the heavy solid curve in Fig. 
n The constraint AiV = 35 to the number of e-folds be- 
tween CMB/LSS and BBO/DECIGO scales tells us that 
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AN 8 " *' 



P Mpi 
from which it follows that 

P t (k g ) _ V{4>g) _ „ 



35, 



(30) 



(31) 



We thus find that the gravitational-wave amplitude at 
DECIGO/BBO scales is 

n GW h 2 = A GW P t (k g ) = A GW rP s {k c )e- 2eAN 
= 1.08 x i -i3 ee -7o e (AJV/35) 
P s (fc c ) \ ( A GW 



2.45 x 10" 



2.74 x 10" 



(32) 



This expression is maximized for e = 1/(2 AN) ~ 1/70 
at a value fl^^h 2 — 5.68 x 10 -16 . Interestingly enough, 
the IGWB detectability through direct detection is max- 
imized for relatively small e, while the detectability with 
the CMB is maximized at larger e j3|| . Given that CMB 
sensitivities are expected to get to r ~ 0.01 in the rela- 
tively near future with the CMB, and then to r ~ 10~ 4 
with a next-generation satellite experiment, it is unlikely 
that this model would produce a direct detection without 
producing a detectable CMB signal. 

Fig-GI a ) shows the region of the n G wh 2 -n t parameter 
space (at BBO/DECIGO scales) that the n s -r parameter 
space shown in Fig. ^ maps to for power-law inflation. 
The breadth in JIgw^ 2 of the region is due to the 30% 
(at 3ct) uncertainty in P s (k c ). If power-law inflation is 
the correct model of inflation, then the IGWB is directly 
detectable with BBO for r > 10~ 3 and r > 10~ 6 with 
DECIGO. 



B. Chaotic Inflation 

In chaotic inflation, the inflaton potential is, 

v{4>) = v 



Ah 



pi 



(33) 



In this family of models, e(</>) = (a 2 /16ir)(Mpi/(f>) 2 , and 
therefore inflation ends when = cf> cn d = <xMpi/(4y / 7r). 
If there are N c e-folds of inflation between CMB horizon 
exit and the end of inflation, then Eq. (|18|) gives us <f> 2 — 
(M|j/167r)(4aiV c + a 2 ). We also have rj(cj>) = a(a - 
1)(M p1 /0) 2 /(8tt) from which it follows that at CMB/LSS 
scales, 



note that the scalar running a s = —2(1 — n s ) 2 /(2 + a) 
is always within the current observational constraints 
since |1 — n s \ < 0.1. This family of models is thus 
parameterized by two parameters: 47 < N c < 62 and 
a < 47V C /19 - 40/19. Note that each choice of (a,N c ) 
maps onto a point in the n s -r parameter space, so we 
could just as well choose n s and r as our two independent 
parameters. If we choose to do so, then we assign iV c and 
a by N c = (l-e)/(l-n a -2e) and a = 4e/(l-n s -2e), 
where e = r/16. 

For a fixed value of N c , this family of models is repre- 
sented by a curve in the n s -r parameter space; a spread 
in the range of values for N c broadens this curve into a 
region in the n s -r parameter space, as indicated by the 
region enclosed by the dotted red curves in Fig. ^ 

Once a and N c are specified, the potential prefactor is 
fixed by 
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(35) 



The gravitational-wave amplitude at direct-detection 
scales is then simply, 



fl G wh 2 



128 V 

T Aow W 



PI 



pi 



(36) 



where the field value 



at the time direct-detection 



scales undergo horizon crossing is given by 



Mpi 
16?r 



(4aiV g + a 2 



(37) 



and where N g = N c — 35 = N c — AN is the number of 
e-folds before the end of inflation that DECIGO/BBO 
scales exit the horizon. For chaotic inflation, the value of 
rit at DECIGO/BBO scales will differ from (and gener- 
ally be larger in amplitude, or more negative than) that 
at CMB/LSS scales. The value of n t at DEICGO/BBO 
scales will differ from that at CMB/LSS scales; it will be 



jiven by n t {4> g 



-2e(<y. 



Ffg-Etb) shows the region of the Vl G -^jh 2 -n t parameter 
space (at BBO/DECIGO scales) that the n s -r parameter 
space shown in Fig. ^ maps to for chaotic inflation. The 
breadth in f^ow^ 2 of the region is due to the spread in 
the a-N c parameter space for fixed P s (k c ) = 2.45 x 10~ 9 ; 
there will be a slight additional vertical broadening be- 
yond that shown due to the uncertainty in this param- 
eter. If chaotic inflation is the correct model of infla- 
tion, then the IGWB is directly detectable with BBO for 
r > 10~ 3 and r > 10~ 6 with DECIGO. 



n s = 1 



4iV c 



(34) 



Noting that 47 < iV c < 62, the constraint n s > 0.9 gives 
us a constraint a < 4iV c /19 — 40/19. The constraint on 
the tensor-to-scalar ratio, r — 16e = 16 a/ (4N C + a) < 1, 
leads to a slightly less stringent limit, a < AN C /15. We 



C. Symmetry Breaking Inflation 

We now consider the Higgs potential, 
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(38) 
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FIG. 3: Regions in the Q.Gwh 2 -nt parameter space for (a) power-law, (b) chaotic, (c) symmetry-breaking, and (d) hybrid 
inflation. The colored (shaded) regions map out the corresponding regions in Fig. Here, the gravitational-wave density 
f^Gwft 2 and spectral index rit are both evaluated at DECIGO/BBO scales. Also shown are the sensitivity goals of BBO and 
DECIGO. 



parameterized by Vo and a Higgs vacuum expectation 
value v. Our treatment of this family of models will par- 
allel that for chaotic inflation. In these models, 4> starts 
near the origin and then rolls toward <f> — u. The slow- 
roll parameters are e(<j>) = (M| 1 2 /4tti/ 4 )[1 - 0/z^) 2 ]~ 2 , 
and n(4>) = (M 2 1 /2ttz/ 2 )[3(0/i/) 2 - 1][1 - (0/^) 2 ]~ 2 , horn 
which we infer an end to inflation, 



fond 



2vr I 



-.1/2 



M 2 



1 + 4tt 



(39) 



The field value at which CMB/LSS scales undergo hori- 
zon crossing during inflation is 



N C M^ 



+ 4% nd - 2u 2 hi(4> end /(j) c ). 



At CMB/LSS scales, 



1-i 
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pi 



1 + 3y c 2 



(40) 



(41) 



where y c = 4> c /v. Since n s is a decreasing function of 
y, the constraint n s > 0.9 requires v > 1.8 Mj>\. The 
prefactor Vo is then fixed by the constraint, 

Vo = ^P s (fc c )(M P1 /^ 2 (i (42) 



V, 
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Once this normalization is fixed, these models are param- 
eterized by v and iV C) and n 8 and r are fixed once these 
two parameters are specified. As in chaotic inflation, we 
may alternatively take as our two free parameters n s and 
r, and then determine v and iV c , although the inversion 
is not as tractable algebraically as in chaotic inflation. 

The gravitational-wave amplitude at direct-detection 
scales is then simply, 



GW 
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— Aow w 
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pi 
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(43) 



where y g 



-jjv, and the field value <p g at the time 



direct-detection scales undergo horizon crossing is given 
by 



^ond 



2 J , 2 ln(</> end /<^), (44) 



and where again, N g — N c ~ 35 = N c — AN is the number 
of e-folds before the end of inflation and the time when 
DECIGO/BBO scales exit the horizon. The value of th at 
DECIGO/BBO scales will differ from that at CMB/LSS 
scales; it will be given by n t (4> g ) = — 2e(0 9 ). We also note 
that the running of the scalar spectral index at CMB/LSS 
scales is 



1 / MpA 2 5 + 3y 2 



V 



(1 



(45) 



We check in our numerical results that all of the models 
we consider are consistent with the bound to this param- 
eter, \a s \ < 0.04. In particular we find that \a s \ < 10~ 3 . 

Fig. Etc) shows the region of the VtQ-^fh 2 -n t parame- 
ter space (at BBO/DECIGO scales) that the n s -r pa- 
rameter space shown in Fig. ^ maps to for symmetry- 
breaking inflation. The breadth in f^ow^ 2 of the region 
is due to the spread in the v—N c parameter space for 
fixed P s (k c ) — 2.45 x 10™ 9 . If symmetry-breaking infla- 
tion is the correct model of inflation, then the IGWB 
will be detectable with BBO and DECIGO. Inciden- 
tally, we have also investigated potentials of the form 
V((j)) = V Q [1 - {(t)/v) p ] 2 with p > 2 [p. In these models, 
the symmetry-breaking scale can be reduced below Mpj, 
however the IGWB amplitude is then reduced below the 
level accessible to BBO/DECIGO for v < 0.1M P i. 



D. Hybrid Inflation 

Hybrid inflation generally requires two scalar fields 
but the phenomenology can be modeled by a sin- 
gle scalar field with the potential, 



V{</>) = V 



(46) 



the potential is not to be taken to be generic within the 
class of hybrid inflation but only as a particular example. 
Other forms exist, such as those found in, e.g., Refs. 
13 EH- Defining y = (f)/fi, we find that the slow-roll pa- 
rameters are given by, e(y) — (Mpi/47r^i 2 )?/ 2 [l + y 2 ]~ 2 , 
and r](y) = (Mpi/47r/i 2 )[l + y 2 ]^ 1 . In these models, e 
is maximized at y = 1 with a value less than unity if 
fi > Mpi/(4 v /7r). Fo r sm aller values of //, inflation ends 
soon after y = 1 [8(| |105| . The dynamics of these mod- 
els resembles those of chaotic inflation, which we have 
already considered, and so we consider them no further. 
New inflationary dynamics arises for \i > Mj>\/ 
and y < 1, and so we restrict our attention here to this 
regime. 

The field value at which CMB /LSS scales undergo hori- 
zon crossing during inflation is, 



Vc. 
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\Vend 
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Vend 
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pi 
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Since y c is taken to be a free parameter, Eq. I|47|l. along 
with 47 < N c < 62, determines the value of y cn d- From 
the slow-roll parameters, 



2y c 2 



2n f i 2 (l + y 2 ) 2 



(48) 



at CMB/LSS scales. The above expression for n s is 
maximized at y = 0, and at this field value becomes 
n s = l+AIp l /2TT/j, 2 , which shows that we can have n s > 1 
in hybrid inflation. The pre-factor Vq is then fixed by the 
constraint, 
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(49) 



Once this normalization is fixed, these models are param- 
eterized by /i and y c , and n s and r are fixed once these 
two parameters are specified. As in chaotic inflation, we 
may alternatively take as our two free parameters n s and 
r, and then determine [i and y c . In particular, these can 
be determined from 



8(n 



1/2 



L) + 2r_ 

2 4(n s -l)+r 
7T [8(n s - l) + 3r] 



1/2 



Mpi. 



(50) 
(51) 



The gravitational-wave amplitude at direct-detection 
scales is then simply given by, 

128 V 
n GW h 2 = — A GWjjr(l + y 2 ), (52) 

where the field value y g at the time direct-detection scales 
undergo horizon crossing is given by 



and the selection of a value for 4> cn d (with the only re- 
quirement that </> C nd > 0). We note that this form for 



y 2 = 21n( 

Vend 



Vend r> 
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(53) 
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and where again, N g = N c — 35 = N c — AN is the number 
of e-folds before the end of inflation and the time when 
DECIGO /BBO scales exit the horizon. The value of nt at 
DECIGO/BBO scales will differ from that at CMB/LSS 
scales; it will be given by n-t(y g ) = —2e(y g ). The running 
of the tensor spectral index, 



1 

4^2 



Ah 



pi 



y 2 (i-y 2 ) 



(i + y 2 ) 4 



(54) 



can be positive in this class of models. Thus, for y < 1, 
at > 0, and the running is positive, indicating that as y 
evolves, the tensor spectral index becomes less negative. 
As we have seen in the previous models, a non-negligible 
gravitational- wave amplitude at CMB/LSS scales leads 
to a "large" amplitude at direct-detection scales pri- 
marily due to a small, negative, tensor spectral index. 
We therefore expect this model to produce the largest 
gravitational- wave amplitude at direct-detection scales. 
We also note that the running of the scalar spectral in- 
dex at CMB/LSS scales is 



2ttV (1 + y, 



2 VI 



(55) 



With the restriction that y c < 1, a s is maximized at 
y = \Jl-°J?>. We note from this that if fi > 0.69M P i 
the observational bound on a s is satisfied for all y c < 1. 
For [i not satisfying this restriction, there will be some 
range of y c which are incompatible with observations. 
This restriction is taken into account in our numerical 
calculations. 

Fig. [3d) shows the region of the Q.G^jh 2 -nt parameter 
space (at BBO/DECIGO scales) that the n s -r parameter 
space shown in Fig. ^ maps to for symmetry-breaking 
inflation. The breadth in flQ-^/h 2 of the region is due to 
the spread in the [J,-N c parameter space for fixed P s (k c ) = 
2.45 x 10~ 9 . If hybrid inflation is the correct model of 
inflation, then the IGWB may be detectable with BBO 
and DECIGO, but it is not guaranteed. 



THE (RUNNING) POWER-LAW 
APPROXIMATION 



During inflation, the value cf> of the scalar field can 
change very little as the scale factor a(t) grows extremely 
rapidly. It is therefore a feature of inflation that a vast 
range of distance scales can correspond to a small change 
in ip. This motivates the power-law expansions (with a 
slight running of the spectral index) in Eqs. (|10|> and 
(|llfl . which assume that the inflaton potential can be ac- 
curately approximated by its Taylor expansion (to second 
order) about a given inflaton value. These power-law ex- 
pansions are particularly appropriate when studyin g th e 
CMB and large-scale structure (e.g., Refs. |35L l4lj l50| ) . 
which involve a spread in distance scales of maybe three 
orders of magnitude. 



However, BBO/DECIGO frequencies are separated 
from those probed by the CMB/LSS by roughly sixteen 
orders of magnitude. The inflaton may thus traverse a 
significant distance, and so it is not obvious that the Tay- 
lor expansion approximation that underlies the power- 
law approximation (even with the running of the spec- 
tral index) will remain valid. For example, in Eqs. i|12|) 
and (|13fl . the tensor and scalar tilt are written in terms 
of the first-order slow-roll parameters, while second- and 
higher-order corrections (e.g., Ref. [83) may be impor- 
tant when extending the power spectrum over large phys- 
ical scales. Similarly, one must also account for higher- 
order derivatives of the tilt, beyond the running consid- 
ered with a s and at- For the calculation performed here, 
higher order corrections are not important as f^cw^ 2 
was directly determined with model parameters describ- 
ing the inflaton potential, rather than through the power 
spectrum. Assuming the Taylor approximation is valid, 
then measurements of P s , n s , and r at CMB/LSS scales 
fix the parameters Pt{k c ), n t , and at in Eq. I|ll|) . which 
can then be used to predict Pt(k g ) oc S1gw^ 2 j the IGWB 
amplitude at BBO/DECIGO scales. An approach based 
on the Taylor expansion was considered in Ref. |35| to 
estimate the GW amplitude at frequencies correspond- 
ing to direct detections. Fig. 0] plots the exact IGWB 
amplitude obtained from the calculation in the previous 
Section versus that obtained from the power-law approx- 
imation. For small IGWB amplitudes, r = 16e — > 0, and 
so the potentials are very close to flat and the power law 
tends to be a good approximation, and is indeed a good 
approximation for the four classes of models we have con- 
sidered. For power-law inflation, where the power spec- 
tra are precisely power laws, the two results are identical. 
For chaotic and symmetry-breaking inflation, e becomes 
large when the IGWB amplitude becomes large, and e 
evolves during inflation in such a way that the power-law 
approximation overestimates the true IGWB amplitude 
at direct-detection scales. The behavior of hybrid in- 
flation is a bit more subtle. The running of the tensor 
spectral index is a t — (r/8)(n s — 1 + r/8) which, for 
n s > 1 — r/8 can be positive. The (running) power-law 
would then suggest that n t will become positive at some 
small distance scale, which cannot be [see Eq. (|13fl]. The 
power-law approximation can then overestimate the true 
IGWB amplitude at BBO/DECIGO scales. On the other 
hand, in hybrid inflation, unlike chaotic and symmetry- 
breaking inflation, e can indeed decrease as inflation pro- 
ceeds, and so the direct-detection amplitude may also be 
underestimated by the power-law approximation. Both 
behaviors are seen in Fig. 21 

It should also be kept in mind that the tensor spectral 
index n t is most generally different at CMB scales than it 
is at direct-detection scales, and it is conceivably measur- 
able at both. Determination of n t at both distance scales 
could therefore distinguish between inflationary models. 
For example, in power-law inflation, nt remains precisely 
constant, while it can change by roughly a factor of two 
for chaotic inflation for models with a directly detectable 
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FIG. 4: Here we plot the inflationary gravitational-wave background amplitude QgwIi 2 obtained with the exact inflationary 
dynamics described in Section IV versus the amplitudes obtained with the power-law approximations (with slow-roll parameters 
fixed by CMB/LSS observations) given in Section lib. The panels show results for (a) power-law, (b) chaotic, (c) symmetry- 
breaking, and (d) hybrid inflation. The regions are models taken from the n s -r parameter space in Fig. (and shown in Fig. 
[iU and the blue dashed curves indicate equality. 



IGWB. Realistically, though, the tensor spectral indices VI. BROKEN SCALE INVARIANT SPECTRUM 

are generically (although not in full generality) small, and 
so running of the tensor spectral index will be difficult to 
measure. 



Finally, the four classes of models we have considered 
are not at all exhaustive, and another inflaton potential 
could yield a direct-detection IGWB amplitude different 
from those we have considered here and different from 
what extrapolation from CMB/LSS would suggest from 
the power-law approximation. For example, in models 
with broken scale invariance [88ll89| . the direct-detection 
amplitude could be considerably different. 



To demonstrate that direct observations of the IGWB 
can distinguish between different forms of the inflaton 
potential, we consider as a toy model the broken scale 
invariant (BSI) potential, which features a sharp change 
in the slope of the inflaton potential at some transition 
scale [32, E^, . Such models have been invoked to ex- 
plain, e.g., the paucity of dwarf galaxies observed around 
the Milky Way H |9ll, HI . 
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Consider a potential of the form, 



V{cj>) = V x 




Act>),<p>0 
cA0),0<O, 



(56) 



where Vq is the overall normalization, A is the slope of the 
potential at CMB/LSS scales and c parameterizes how 
the slope changes after the break at <fi = 0. We allow 4> c 
to be a free parameter, only requiring that it be before the 
break in the potential at <f> = 0. This freedom supposes 
that field value at which inflation ends is not necessarily 
determined by the form of the potential in Eq. (|56|l . In 
order to choose 4> c , we place the break (i.e., (f> = 0) N 
e-folds from cj> c . A natural choice for No is 10, since 
we suppose that CMB/LSS scales constrain the inflaton 
potential from 10 4 Mpc to 1 Mpc. The normalization of 
the scalar power-spectrum then fixes the normalization 
Vq of the potential through the expression, 
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We then integrate Eq. I|18[l. assuming the transition at 
4> = has a negligible contribution, between <p c and <p g 
with N = 35 in order to find <f> g . We require that inflation 
not end before we reach <fr g . For V(4>) as in Eq. (|56H . 
we find that inflation ends soon after the field reaches 
a value, 0* = Mpi/(4y / 7r) — (cA)^ 1 . This then places a 
constraint on the combination cA, 



cA < 



4^ 



Mpi VI + 4(iV ff - JV ) 
At CMB/LSS scales we find that r)(<f) c ) = 0, and 



<<t>c) 



A 2 Mp i 



AA 2 M 2 y N + 16tt' 



(58) 



(59) 



From the above expression and Eqs. (|12fl and (|16fl . we can 
see explicitly that n s and r only depend on our choice of 
N and A. At BBO/DECIGO scales we find 



9 ,^9 9 JA 2 c 2 M 2 (N Q - N a ) + 4tt 

where N g is the number of e- foldings between </> c and <j) g . 
In order for there to be N g e-folds between <fr c and 4> g the 
slope of the potential cannot be too large, requiring 



cA < 



M vl ^N g - N 



(61) 



Comparing this to Eq. I|58|) . we find that this constraint 
is slightly less restrictive. We can see that this amplitude 
depends not only on N and A but also on c. Therefore, 
potentials that share approximately the same Taylor ex- 
pansion at CMB/LSS scales, but different expansions at 
BBO/DECIGO scales, will produce overlapping observa- 
tions in the (n s , r) plane at CMB/LSS scales and different 



gravitational-wave amplitudes at BBO/DECIGO scales. 
With the constraint in Eq. (J5BJ, we find that as c in- 
creases towards its maximum value (for a fixed A), the 
amplitude of the IGWB changes by an order of magni- 
tude. For example, for n s = 0.9 and r = 0.27 we have 
1.0 x 10~ 16 < fl GW h 2 < 1.0 x 10~ 15 ; for n s = 0.99 
and r = 3.16 x 10~ 3 we have 2.1 x 10~ 18 < ^ G w^ 2 < 



2.1 x 10 17 ; and for n s — 1.0 and r = 
have 2.1 x 10~ 20 < n GW h 2 < 2.1 x 10" 



= 3.18 x 10~ 5 we 
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VII. DISCUSSION 

In this paper we have calculated the gravitational- 
wave amplitude at direct-detection scales for four classes 
of inflationary potentials with parameters consistent 
with current constraints from the CMB and LSS. The 
gravitational-wave amplitude f^cw^ 2 is proportional to 
the height V(<fi g ) of the inflaton potential at the time that 
direct-detection comoving scales exit the horizon. Our 
current theoretical understanding does not fix V((j) g ); it 
is constrained to be V 1/4 < 3.4 x 10 16 GeV from the 
CMB, and it could conceivably be as low as T ~ 1 MeV 
without violating observational constraints. Moreover, 
detectability of the IGWB with BBO or DECIGO re- 
quires 

yl/4 > 1Q 15 Q eVj dose tQ the 

upper allowed limit . 
It thus seems, a priori, that detectable models occupy a 
small region of parameter space. 

That said, however, there are indeed constraints 
to inflationary models that come from the CMB and 
large-scale structure, notably constraints to the density- 
perturbation amplitude and spectral index. Fig. indi- 
cates that when we go through the exercise of writing 
down simple functional forms for the inflationary poten- 
tials and imposing current constraints, there are large re- 
gions of parameter space that lead to directly detectable 
IGWB amplitudes. In particular, for the symmetry- 
breaking potential, which looks perhaps like the type of 
Higgs potentials we might associated with grand unifi- 
cation, current constraints lead to a directly detectable 
IGWB amplitude. 

The promise of detectability traces back to the fact 
that SIqw'* 2 oc V oc (V) 4 ^ 3 , the last proportional- 
ity tracing back to Eqs. © and © for fixed density- 
perturbation amplitude P s . Thus, if the potential is 
extremely flat, V -> 0, then the IGWB will be tiny. 
However, if the potential takes a form in which V' ~ 
V/<p, which seems theoretically natural, then the re- 
quired density-perturbation amplitude is achieved with 
V ~ i()( 15 ~ 16 ) GeV, the range that produces an accessi- 
ble IGWB amplitude. 

There is of course still plenty of room for inflation to 
be correct and for the IGWB amplitude to be well below 
the BBO or DECIGO threshold. For example, in power- 
law inflation and chaotic inflation, the IGWB amplitude 
becomes small when n s — > 1; i.e., when scale invariance 
is achieved which, in these models corresponds to small 
V' . On the other hand, a value n s — > 1 does not, more 
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generally, imply a small IGWB amplitude. For example, 
in hybrid inflation one can have n s = 1 if = 1/2 [cf., 
Eq. Ij48(l ]. and for [i > 1.8 Mpi, the potential can reach 
values at CMB/LSS scales of V ~ 3 x 10 16 GcV, which 
even after the decrease to BBO/DECIGO scales remains 
within reach of detection, as shown in Fig, ffid). 

There may of course be alternatives to inflation, 
such as cyclic models |93j or the pre big-bang model 
[H HI HE US HI, that have completely different 
IGWB spectra. Although the cyclic model predicts a 
blue tensor spectrum, which might improve detectabil- 
ity at small scales, BBN constrains the amplitude of the 
gravitational-wave amplitude to be or ders o f magnitude 
below the BBO/DECIGO sensitivities [Tonj . 

Our conclusion is that direct detection of the IGWB 
is unlikely without detection with the CMB polarization. 
Still, direct detection could be extraordinarily valuable 
even if the IGWB is detected first in the CMB. Direct 
detection would provide yet another cross-check that the 
curl component in the CMB polarization is due to grav- 
itational waves, as opposed to some other mechanism 
(e.g., vector modes, cosmic shear, or foregrounds). Since 
it occurs on such vastly different distance scales, direct 
detection can verify that it is a nearly scale-invariant 
spectrum of gravitational waves, as predicted by infla- 
tion, as opposed to some other phenomenon that might 
only produce large- wavelength gravitational waves. It 
would provide evidence for the continuation of infla- 
tion to distance scales well beyond those implied by the 
smoothness of the Universe suggested by the successes 
of BBN. The large lever arm provides an opportunity to 



discriminate between inflationary models that produce 
the same CMB/LSS observables. Even within the con- 
text of a given inflationary potential, the large lever arm 
associated with direct detection may allow a measure- 
ment of inflationary parameters that may be more pre- 
cise than those accessible with the CMB/LSS alone. For 
example, an uncertainty of 10% in e from the CMB/LSS 
translates to a ~ (10 15 ) 2 ~ 1000 uncertainty in the 
BBO/DECIGO amplitude. Thus, a detection alone, with 
no better than an order-unity amplitude uncertainty, cor- 
responds to a measurement of e to roughly 0.02, probably 
better than is accessible with the CMB/LSS alone. Fi- 
nally, the deci-Hertz IGWB amplitude counts the num- 
ber of relativistic degrees of freedom at temperatures a 
bit above the electroweak symmetry-breaking scale, and 
may thus be used to probe for new degrees of freedom as- 
sociated with supersym metry of some other new physics 
at the electroweak scale |101| . The direct search for in- 
flationary gravitational waves may thus be warranted. 
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